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Variable structure control system design methods are applied to the problem of aircraft 
spin recovery. A variable structure control law typically has two phases of operation. The 
reaching mode phase uses a nonlinear relay control strategy to drive the system trajectory to 
a pre-defined switching surface within the motion state space. The sliding mode phase 
involves motion along the surface as the system moves toward an equilibrium or critical 
point. Analysis results presented in this paper reveal that the conventional method for spin 
recovery can be interpreted as a variable structure controller with a switching surface 
defined at zero yaw rate. Application of Lyapunov stability methods show that deflecting 
the ailerons in the direction of the spin helps to insure that this switching surface is stable. 
Flight test results, obtained using an instrumented hand-launched glider, are used to verify 
stability of the reaching mode dynamics. 
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Nomenclature 

wing span 

mean aerodynamic chord 
angular momentum vector 
roll moment of inertia 
pitch moment of inertia 
yaw moment of inertia 
relay gain 

total rolling moment 

aerodynamic rolling moment 

rolling moment due to aileron deflection 

total pitching moment 

aerodynamic pitching moment 

pitching moment due to elevator deflection 

longitudinal magnetic field strength 

lateral magnetic field strength 

vertical magnetic field strength 

total yawing moment 

aerodynamic yawing moment 

yawing moment due to rudder deflection 

roll rate 

pitch rate 

dynamic pressure 

yaw rate 

servo rate limit 

reference area 
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t = time 

u = control vector 

V = Lyapunov function 

x = state vector 

y = subspace vector 

A = control law update rate 

8 4 = aileron deflection 

S E = elevator deflection 

S R = rudder deflection 

£ = servo deadband 

a = switching surface 

V = control deflection noise 

I. Introduction 

I T has long been known that the most effective spin recovery technique for fuselage-heavy aircraft is to use the 
ailerons to roll the aircraft into the spin. This technique was first proven effective by A.I. Neihouse in 1940, using 
scale models in the free-spinning wind tunnel at Langley Aeronautical Laboratory. 1 This recovery strategy might be 
considered counterintuitive because the pro-spin aileron deflection will tend to aggravate the roll rate component of 
the aircraft's angular momentum. The objective of this paper is to explore this well-known spin recovery strategy 
using a recently-developed integrated toolset for modeling and analysis of highly dynamic aircraft motions 

The integrated toolset combines existing nonlinear methods for aerodynamic modeling, stability analysis, and 
control system design. The nonlinear lifting line aerodynamic modeling method was first introduced by Piszkin and 
Levinsky in 1975. 2 * This method differs from conventional lifting-line methods in that it allows the section lift 
coefficient to be expressed as a nonlinear function of angle-of-attack. This extension yields accurate aerodynamic 
information over a wide angle-of-attack range. 

Hsu introduced the cell-to-cell mapping method for nonlinear system stability analysis in 1987/ The power of 
the cell-to-cell mapping theory comes from the fact that the map itself represents a linear Markov chain. 
Consequently, linear system stability analysis tools can be applied to the map and the results reveal how the 
underlying nonlinear system will behave. 

Variable Structure Control (VSC) is a high-speed switching feedback control technique that has been developed 
for control of nonlinear dynamic systems. 4 This method was added to the integrated toolset because the relay 
control behavior of VSC is very similar to automatic upset recovery systems used on full-scale aircraft. For 
example, some fighter aircraft are equipped with a spin arrow display that requires the pilot to apply full control 
authority in the direction of the spin arrow. 

The nonlinear lifting line aerodynamic modeling method and the subspace cell-to-cell mapping method are 
described in greater detail in a recent contractor report/ The contractor report also provides a user guide to the 
software that has been developed for each method. These software programs were used to generate the results 
described in this paper. 

A brief review of the subspace cell-to-cell mapping method is provided in Section II of this paper. Section III 
describes how coventional spin recovery control strategies are related to variable structure control laws. Section IV 
describes the hand-launched glider that was used to gather flight test data. Section V presents analysis and flight- 
test results. Several important conclusions are offered in Section VI. 

II. Subspace Cell-to-Cell Mapping 

The cell-to-cell map represents the transition probability of solution trajectories moving from one cell to another. 
The method works well for systems of dynamic order less than about five. When the dynamic order exceeds five, 
the number of cells needed to accurately represent the system dynamics becomes prohibitively large. The approach 
used in this research has been to define a two-dimensional subspace within the state-space domain of interest. The 
cell-to-cell mapping method is then used to create a cell map on the two-dimensional subspace. 

The subspace mapping approach is illustrated in Fig. 1 for a third-order system. The state variables are Xi, x 2 , 
and x 3 . The two-dimensional subspace is defined by the variables y 3 and y 2 , such that Xj = yi, x 2 = y 2 , and x 3 = 0. 
The two-dimensional subspace is divided into twenty-four cells in Fig. 1. The cell-to-cell map for this system is 
constructed by starting several trajectories from within each cell. These trajectories are propagated forward for a 
fixed time period. The point at the end of the trajectory is projected onto the two-dimensional subspace and the cell 
in which this point lands is recorded. As Fig. 1 illustrates, solutions that start from the one cell do not necessarily 
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land in the same cell. These differences are 
captured in the transition probabilities that 
form the cell-to-cell map. For example, if 
100 trajectories are started from cell #2, and 
40 of those ends in cell #4, then the (4,2) 
element of the cell-to-cell map is set to 
40/100. 

The cell-to-cell map that is created on the 
two-dimensional subspace is obviously not 
sufficient to capture the stability properties of 
the entire state-space. However, significant 
insight can be gained by locating and 
orienting the subspace near features of 
interest. One example is to study the 
behavior of the system near a critical point by 
defining several orthogonal subspaces that 
intersect the critical point. 

In analyzing cell-to-cell map properties, Hsu devised ways to sort and categorize how solution trajectories enter 
and leave the cells. Two of the categories that are most relevant to this research are absorbing cells and transient 
cells. A transient cell has the property that once a solution leaves the cell, it does not return. An absorbing cell has 
the property that once a solution enters the cell, it does not leave. Stable critical points in the nonlinear system state- 
space are typically found near one or more absorbing cells. 

It is possible to locate the absorbing cells by studying the long-term behavior of the Markov chain. The list of 
transient cells that have a high probability of moving into an absorbing cell determine the basin of attraction for the 
absorbing cell. Solutions that start in transient cells located near the edge of an attraction basin would presumably 
take a longer time to reach the absorbing cells. The expected absorption time is used to quantify this relationship. 
Absorbing and transient cells will be identified in subsequent analysis results using plots of the expected absorption 
time from within the two-dimensional subspace. 

III. Variable Structure Spin Control 

The VSC control strategy is divided into two phases. The first phase of the control strategy is designed to drive 
the state trajectory onto a specified subspace called the switching surface. Behavior of the closed-loop system 
during this phase is often called the reaching mode. The second phase of the VSC control strategy defines the 
system behavior while the system is operating on or near the switching surface. This phase is called the sliding 
mode. A linear state feedback structure is generally used for sliding mode control. Its objective is to drive the 
system state to a desired equilibrium point by "sliding" the trajectory along the switching surface. 

This paper will focus on the reaching mode phase of the VSC control technique. A constant gain relay is one of 
the most common elements used to implement reaching mode control laws. Assume that the function o(x) = 0 
represents a switching surface. If x is the state vector and u is the control variable, the constant gain relay control 
law is given in Eqs. (1) and (2). 



/ 

Xi =Yi 

Figure 1. Subspace Cell-to-Cell Map Construction 


u = ksgn [(J(x) ] 


(1) 


sgn[ a(x)] 


(+1, <7(x)> 0 
{-1, a(x)<0 


(2) 


The design philosophy of the VSC reaching mode is to insure that the domain of attraction to the switching 
surface is as large as possible. This objective is achieved by choosing switching surfaces such that asymptotic 
stability can be guaranteed using the Lyapunov's direct method. A Lyapunov function V(x) is defined in terms of the 
switching surface, as shown in Eq. (3). The Lyapunov function V(x) is positive definite for all nonzero x and V(0) = 
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0 is satisfied as long as (7(0) = 0. Under these conditions, the system state trajectory will converge to the switching 
surface as long as the control law is designed to keep dV/dt < 0. The time derivative of V(x) yields the global 
stability condition relative to the switching surface. 

The reaching mode Lyapunov stability condition, expressed in Eq. (4), will be used in this paper to investigate 
control strategies for aircraft spin recovery. The objective of the VSC control law will be to drive the aircraft away 
from the spin critical point and toward a stable switching surface. For aircraft spin control, it would seem logical to 
design switching surfaces that are related to the aircraft's angular momentum. Relative to a conventional body-fixed 
axis system, the aircraft's angular momentum vector is written as shown in Eq. (5). 6 Note that it has been assumed 
that the cross-products of inertia are small. 

Separate switching surfaces are defined for each component of the angular momentum vector in Eqs. (6) - (8). 
The associated Lyapunov stability condition for each switching surface is listed in Eqs. (9) - (11). Note that cross- 
axis control effects, such as yawing moment due to aileron deflection, have been ignored for simplicity. 


V(x) = <r(x)/2 (3) 

dV/dt = (7da/dt < 0 (4) 

h=[I x p ly q I : r] T (5) 

<7,=I x p (6) 

<72 = ly q ( 7 ) 

(Ts=I z r (8) 

dVj/dt = a / d<7i/dt = I x p [ (I y - 1-) qr + q har S b (L A + L SA S A ) ] (9) 


dV/dt = o 2 daVdt = I y q [ (I. - Ij pr + q bar S c bar (M A + M m 8,) ] 


( 10 ) 


dV 3 /dt = a, da 3 /dt =I : r [ (I x - I y ) pq + q bar S b (N A + N S r Sr) ] 


( 11 ) 


Spin recovery strategies can be developed by considering the effect of constant gain relay control laws on the 
stability conditions defined above. First consider Eq. (12), which is a constant gain relay control law for the rudder. 
The gain k 3 is set equal to the maximum rudder deflection S Rmax . This control strategy defines a switching surface <7 3 
= 0 and is governed by the Lyapunov function V 3 (x). In a right spin, for example, the yaw rate is positive. 
Therefore, (J 3 is positive in a right spin and the rudder control law will set Sr = +S Rmax . A negative yawing moment 
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results because the rudder sign convention is such that N ( ,- K is negative. This control moment helps to force dV3/dt 
negative. If dV^/dt remains negative, the spinning aircraft will be driven to the a 3 = 0 switching surface, which has 
the net result of driving yaw rate to zero. This rudder relay control law mimics the conventional wisdom that the 
rudder should be deflected to oppose the spin. 


S R = k 3 sgn [ (7 3 (x)J 


( 12 ) 


Now consider the constant gain control law for the aileron defined in Eq. (13). This aileron relay control law 
defines a switching surface at < 7 ; = 0. The aileron relay control gain will first be chosen such that ki = - S Amax • A 
sign convention of Lg A > 0 implies that this aileron relay control law will deflect the aileron to generate a rolling 
moment in a direction opposing the spin. For example, both the roll rate p and yaw rate r will be positive in a right 
spin. The aileron relay control law will therefore command S A = -S Amax in a right spin. This command deflection 
will produce a negative rolling moment, which will help to force dV/dt negative. 


S A = k, sgn [ G t (x) ] 


(13) 


It would seem logical that a combination of the rudder and aileron relay control laws described above would 
yield an even more effective spin recovery technique. The rudder control law will tend to drive the airplane to the 
G } = 0 surface and the aileron control law will tend to drive the airplane to the <J> = 0 surface. Experimental results, 
shown later in Section IV, reveal that these two objectives conflict, causing the airplane to oscillate between the 
surfaces rather than actually reach either one. The aircraft eventually reaches the intersection of the surfaces, where 
a, = o> = 0, but the technique does not lead to the most rapid recovery. 

The aileron will be deflected into the spin if the relay control law gain is set such that k t = S Amax . In a right spin, 
with roll rate positive, the aileron relay control law will set S A = S Amax . Because Lg A is positive, this aileron setting 
will produce a positive rolling moment. A positive rolling moment will introduce a positive contribution to dVj/dt. 
This choice in aileron control gain will tend to force the aircraft away from the <7, = 0 switching surface. 

The value in driving the aircraft away from the a ) = 0 switching surface is seen by careful consideration of the 
inertia coupling terms in the Lyapunov stability conditions. The inertial coupling term in the dV 3 /dt stability 
condition is (I x - I v ) pqr. For a fuselage-heavy airplane, I Y is larger than I x , which means that the inertial coupling 
term will be negative, in a right spin, as long as the pitch rate q is positive. Therefore, a beneficial effect is achieved 
with regard to the stability of the a 3 — 0 switching surface if q is forced positive in a right spin. Ignoring the 
aerodynamic pitching moment and assuming neutral elevator, the pitch acceleration equation is given by Eq. (14). 
Forcing the aircraft away from the a ) = 0 switching surface tends to increase the aircraft roll rate. Equation (14) 
shows that the increased roll rate will introduce pitch acceleration in the same direction. Consequently, the pitch 
rate will become positive in a right spin when the ailerons are deflected into the spin. 


/,, dq/dt = (I, - IJ P r 


(14) 


The spin recovery technique discovered by Neihouse is therefore identical to a VSC strategy that centers on the 
a 3 = 0 switching surface. The rudder and aileron are both used to insure stable reaching mode dynamics relative to 
this surface. The rudder relay control law provides a direct negative contribution to the dV 3 /dt stability condition. 
The pro-spin aileron relay control law provides an additional negative contribution indirectly, through inertial 
coupling. As a result, the airplane is expected to rapidly reach the <J 3 = 0 switching surface. After the stable 
switching surface is reached, the sliding mode phase of the VSC strategy would be used to complete the spin 
recovery. 
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IV. Vehicle Description and Modeling 

A small hand-launched glider was constructed to provide test data and verify the theoretical results. The glider, 
shown in Fig. 2, was built by modifying a powered, remote control, aerobatic model airplane kit. The main structure 
of the glider is made from expanded-polypropylene foam reinforced with carbon-fiber wire. A battery tray 
containing four AAA batteries is mounted in the slot originally intended for the remote-control receiver batteries. 
The center-of-gravity for the flight-ready glider is located 20 cm behind the nose of the aircraft and is close to the 
fuselage centerline. The moments of inertia for the glider were estimated by measuring the weight and physical 
dimensions of each component before final construction. The heaviest components are the fuselage (74 g), battery 
tray (56 g), and the printed circuit board (32 g). The estimated mass and inertia characteristics are listed in Table 1. 

A single experimenter can perform all of 
the steps necessary to flight test the glider. 

Setscrews on each servo control arm are Table 1. Glider Geometry and Mass Properties 
used to center each control surface. Once 
power is connected, a green LED 
illuminates to serve as a power indicator. A 
handheld computer can be attached to the 
printed circuit board using a four-pin 
modular connector. A software program 
written for the handheld computer allows 
the experimenter to download and store data 
from the airplane. It can also be used to 
change the control sequence to be flown 
during the next test. 

A flight test begins when the 
experimenter presses a momentary reset switch located on the circuit board. A yellow LED illuminates to indicate 
that the circuit board is now recording data. A pulse input is applied to the left aileron after a two-second delay. 
This aileron pulse signals the experimenter to launch the glider. The on-board microcontroller will begin execution 
of the programmed control sequence two seconds after the aileron pulse. Each flight lasts from five to eight seconds, 
depending on the control sequence, initial velocity, and initial altitude. The yellow LED is automatically turned off, 
and data is no longer collected, after a time period of fifteen seconds. Subsequent test flights can be initiated by 
pressing the reset switch again. Data 
collected from each flight is appended to the 
data from previous flights. More than sixty 
flights can be stored before the handheld 
computer is needed to download the data and 
to reset the flash memory. 

The glider is equipped with an 
application-specific instrumentation board. 

The printed circuit board was developed to 
collect data and control the glider in flight. 

The circuit is controlled by a microcontroller 
that is driven by a 20 MHz crystal and 
includes a multi-channel 10-bit analog-to- 
digital (A/D) converter. The microcontroller 
drives up to four servos. For this application, 
micro-miniature servos are used to control the 
left aileron, right aileron, rudder, and 
elevator. 

A yaw rate gyro is mounted near the 
center of the circuit board. This gyro has a 
dynamic range of +300 deg/s. Using the 10- 
bit A/D converter of the microcontroller 
results in a yaw rate resolution of better than 
1 deg/s. Magnetic field strength 

measurements are obtained using three Figure 2. Hand-Launched Glider 



Wing Span, b 76.5 cm 

Length 65.0 cm 

Wing Reference Area, S 1400 cm 2 

Mean Aerodynamic Chord, c bar 18.5 cm 

Aspect Ratio 4.2 

Mass 0.226 kg 

Roll Moment of Inertia, I x 16.2 kg-cm 2 

Pitch Moment of Inertia, I y 31.8 kg-cnr 

Yaw Moment of Inertia, I z 45.4 kg-cm 2 
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magnetoresistive sensors. Single-supply instrumentation amplifiers amplify the differential voltage output from 
these sensors. This combination leads to a minimum resolution of about 150 nano-Teslas. The combination of three 
magnetoresistive sensors forms a three-axis magnetometer. The magnetometer was calibrated by comparing its 
response to the IGRF2000 geodetic model published by the National Geographic Service. 

The VSC spin recovery control laws require estimates of sgn[o/x)] and sgn[<J 3 (x)J. The sign of the yaw rate 
gyro signal is taken as the value of sgn[o 3 (x)]. The value of sgn[o/x)] is estimated using the magnetoresistive 
sensors. The rate of change in the lateral magnetoresistive signal is related to the body-axis angular rates through 
Eq. (15). Because r is already measured by the yaw rate gyro, Eq. (15) can be solved for roll rate p. The sign of the 
roll rate estimate in Eq. (16) is assigned to sgn[o 1 (x) ]. A backward difference approximation was used to estimate 
the magnetometer signal rate dm/dt. The backward difference expression shown in Eq. (17) where m v (k) is the 
current lateral magnetic sensor reading, m y (k-l) is the previous sensor reading, and A is the time difference between 
readings. 


dm/dt = -r m x + p m z 


(15) 


sgn[p] = sgn [dm/dt + r mj sgn[m~] 


(16) 


dm v (k)/dt = [m/k) - m/k-l)]/A 


(17) 


Sensor data are stored on a flash memory chip. The nonvolatile flash memory chip is organized with 2048 
pages of 264 bytes each (540 KB total). Each data record contains twelve bytes and data records are written at a rate 
of 40 records per second. As a result, sensor and control data are written to the flash memory chip at a rate of 480 
bytes per second. The forward difference approximation in Eq. (17) is updated at a rate of only 4 times per second 
(A = 0.25 s). The VSC control laws are updated at a rate of 4 times per second as well. 

The dynamic model for the control servos is shown in Fig. 3. During bench testing, it was noted that the servos 
responded to maximum amplitude position commands by moving at a nearly constant rate. The constant control rate 
R = 2.8 rad/s (160 deg/s) was set by observing that each control surface travels forty degrees in 0.25 s. A deadband 
of s= 0.03 rad (2 deg) is also assumed. A random noise signal is added to each control surface position for post- 
flight data processing. This zero-mean, Gaussian noise signal v is intended to model nonlinearity in the control 
surface linkage caused by control-rod flexing, friction, and control horn backlash. The intensity of the control 
surface noise is set to 0.01 s" 1 . 

A nonlinear filtering method is used to 
reconstruct motion state variables that are not 
measured in flight. Extended Kalman Filters 
(EKF) have been traditionally used for this 
purpose. The EKF linearizes the aircraft 
motion model and measurement functions, 
relative to the current mean state, so that a 
linear Kalman filter can be utilized. Its 
obvious limitations are that the motion 
dynamics must be predominately linear and 
the noise sources must be Gaussian. These limitations make the EKF less suitable for reconstructing state 
trajectories from highly dynamic maneuvers. To improve performance, the Unscented Kalman filter (UKF) was 
adopted for this application. The UKF utilizes a set of special vectors, called sigma-points, to represent the state 
probability distribution. 7 The UKF is very similar to the square-root implementation of the EKF because the sigma- 
points are constructed using the square root of the state covariance matrix. The computational burden of the UKF is 
therefore similar to the EKF, which makes it an ideal substitute for the EKF. Standard equations of motion are used 
to represent the glider dynamics and quaternions are used to describe the aircraft's attitude. 



Figure 3. Servo Model 
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V. Analysis and Flight Test Results 

The nonlinear lifting line software was used to estimate a complete aerodynamic model for the hand-launched 
glider. The wing was divided into two distinct aerodynamic surfaces and the vortex structure for each wing surface 
is solved separately. Each wing surface is divided into seven sections and five rows of wake vortices are used for 
each surface. The horizontal tail is also split into two separate surfaces, with four sections for each surface. The 
horizontal tail surfaces are affected by the wake of the wing surfaces. The vertical tail has four sections as well. 
The wing airfoil section is the NACA 00 12. 8 The tail surfaces are modeled using the NACA 0006 airfoil section. 

The estimated lift coefficient characteristics predict that stall occurs at approximately twenty-four degrees angle- 
of-attack. The static margin is about 0.12. Conventional flight dynamic modes are readily identified using spectral 
analysis of the glider model linearized about an equilibrium glide condition. Figure 4 shows the glider characteristic 
roots as a function of angle-of-attack for a constant airspeed of 6 m/s. Each dot represents a two-degree increment 
in the range from six to thirty degrees angle-of-attack. 

For smaller angles-of-attack, the roll mode time 
constant is about 0.05 s and the spiral mode time 
constant is 50.0 s. The Dutch-roll mode has a 
natural frequency of 4.6 rad/s and a damping ratio 
of 0.6. The phugoid mode is stable, and has a 
natural frequency of 1.4 rad/s with a damping ratio 
of 0.4. Two stable real roots, located at -2.2 and - 
14.2 rad/s, define the short-period mode. 

The aircraft remains stable as the angle-of- 
attack is increased beyond stall. The characteristic 
roots associated with the two longitudinal modes 
move toward one another. The natural frequency 
of the phugoid mode increases to about 2.5 rad/s 
while the short-period real roots merge to form a 
set of complex-conjugate roots with natural 
frequency of 9.1 rad/s and damping ratio 0.45. The 
Dutch roll mode roots move slightly closer to the 
origin. The time constant of the roll mode is 
increased and the time constant of the spiral mode 
root is decreased. 

Data from over one hundred test flights of the 
hand-launched glider was collected. The test 
flights were conducted over a period of several 
weeks, in calm winds, and at temperatures above 
70°F. Standard sea level conditions were assumed 
for all data processing because the test field elevation was approximately 60 m above sea level. Data collected 
during each flight was post-processed using the UKF method described previously. Solutions were obtained for 
approximately one-half of the test flights. The majority of the unsuccessful cases resulted when the state covariance 
matrix of the UKF became negative definite. This problem generally occurred at very large sideslip angles and 
might be caused by inadequate aerodynamic modeling of the glider fuselage. 

Flight testing was conducted in two phases. The first phase involved a series of functional test flights to verify 
sensor calibration and control surface movement. These tests were also used to compare the recorded sensor 
readings to predictions from the aerodynamic model. The control deflection inputs for the functional test flights 
included aileron doublets, rudder doublets, and elevator pulses. The second phase of flight testing was devoted to 
the evaluation of the VSC reaching mode control laws. 

Figure 5 shows an example functional test flight result. The glider response to an aileron doublet input is 
depicted in this figure. The solid line is the UKF prediction while the dots represent the sensor measurements. The 
dashed line is the aileron command provided by the on-board microcontroller. The aileron doublet is applied four 
seconds into the flight. The aircraft immediately rolls to a nearly inverted attitude as a result of the doublet input. 
The vertical magnetic and yaw rate sensor readings are shown as dots in the lower half of Fig. 5. These 
measurements are reasonably close to the response predicted by the aerodynamic model and the UKF. Angular rates 
on the order of 100 deg/s are measured by the yaw rate gyro and predicted by the model. The response of the 
vertically mounted magnetic sensor is also shown to be comparable with model predictions. 



Figure 4. Glider Characteristic Roots 
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Figure 5. Aileron Doublet Response 


The subspace cell-to-cell mapping analysis of the spin recovery control strategies will focus on a two- 
dimensional subspace that is aligned with yaw rate and roll rate. The absorption times for cells in this subspace are 
shown in Fig. 6a for neutral control settings and Fig. 6b for spin entry controls. The transient cells are shaded in 
Figure 6 such that a darker shade of blue represents a longer absorption time. Two absorbing cells are found near 
the origin of Fig. 6a. They are left unshaded in Fig. 6a. These absorbing cells are associated with the constant speed 
equilibrium glide condition. The spin entry assumes control settings for a right spin: elevator 20 degrees trailing 
edge up, ailerons 20 degrees right wing down, and rudder 20 degrees trailing edge right. A single absorbing cell is 
found in Fig. 6b. This absorbing cell locates the predicted conditions for a steady right spin. 

The short flight duration and low initial altitude prevented the aircraft from settling in a clearly defined steady 
spin condition. However, flight test results confirm that the glider was able to reach the critical spin conditions 
predicted by the cell maps. Two spin entry trajectories, obtained from glider flight testing, are shown in Fig. 7. The 
open circles in these plots mark where spin entry controls are applied. The open squares denote where the spin 
recovery control laws are initiated. Both trajectories are clearly attracted to the location of the steady spin absorbing 
cells. The path of the two trajectories appears, at first glance, to be vastly different. Flowever, closer inspection 
reveals very similar characteristics in both trajectories. The differences in path are likely due to slightly different 
initial conditions when the spin entry controls are first applied. 
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a) Neutral Controls 


b) Spin Entry Controls 


Figure 6. Cell Maps for Neutral and Spin Entry Controls 



Yaw Rate, rad/s 


Figure 7. Spin Entry Trajectories from Flight Test 


Yaw Rate, rad/s 


Two VSC recovery control laws are evaluated in the next series of figures. The first recovery strategy uses a 
constant gain relay rudder control to oppose the spin and a constant gain relay aileron control, also deflected to 
oppose the spin. This strategy, denoted RAA (rudder+aileron against), defines control deflections in the direction of 
the switching surfaces <7/ = 0 and <7j = 0. The second strategy, denoted RAW (rudder+aileron with), has the aileron 
deflected in the direction of the spin, or against the switching surface G ) = 0. The second strategy mimics the 
conventional method of spin recovery. The elevator is returned to the neutral position in both cases. 

Figure 8 presents example flight test results for the RAA spin recovery control law along with the associated cell 
absorption probability plot. The open square in Fig. 8b marks the point at which the spin recovery control law is 
activated. The open triangle marks the end of the flight. The cell-to-cell map analysis reveals two unshaded 
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absorbing cells. They are designated by the unshaded cells in Fig. 8a. The transient cells in Fig. 8a are shaded so 
that the darkest blue shades denote 2 steps of 0.35 second duration. As a result, the glider is expected to recover 
from the spin in less than one second. The flight test result confirms this result, but also shows that the path to 
recovery is not necessarily direct. 

Flight test results for the RAW spin recovery control strategy are presented in Fig. 9. The plot of cell absorption 
times for this control law reveals two unshaded absorbing cells. These cells are located on the <T ; = 0 switching 
surface with roll rates of approximately 3.5 rad/s (200 deg/s). The transient cells shown in Fig. 9a are shaded such 
that the darkest shade represents 3 steps of 0.12 second duration. Consequently, the glider is expected to move from 
the spin critical point to the nearest absorbing cell in less than one-half second. Figure 9b demonstrates that the 
glider response trajectory moves to the positive roll rate absorbing cell. The trajectory remains near the absorbing 
cell for a short period of time, before returning to the origin as the glider reaches the ground. 

The effectiveness of the RAW spin recovery strategy is apparent in the response trajectories presented in Fig. 10. 
For comparison, time histories from the RAA flight test are shown in Fig. 10a and time histories from the RAW test 
are shown in Fig. 10b. The open circles in these plots represent the time at which the spin entry controls are applied. 
The open squares mark the time at which the spin recovery control is activated. Figure 10b shows that the pitch rate 
remains positive once the RAW recovery is applied. This characteristic was shown in earlier analysis to help the 
airplane move toward the <J; = 0 switching surface. The RAW control law forces the aircraft into a steady roll. 
Figure 9b demonstrates that a steady roll rate of about 200 deg/s is achieved in less than one-half second from the 
time that the RAW control is activated. The RAA control law does not appear to effectively reduce angle-of-attack. 
in fact, the glider angle-of-attack exceeds twenty degrees for nearly one second after the RAA control is activated. 
The angle-of-attack is fairly small when the RAW control law is activated, but one can see that this control law 
keeps the angle-of-attack close to zero throughout the remainder of the flight. 



Yaw Rate, rad/s 


Yaw Rate, rad/s 


a) Cell Map 


b) Flight Test Trajectory 


Figure 8. RAA Spin Recovery Cell Map and Flight Test Trajectory 
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a) Cell Map b) Flight Test Trajectory 

Figure 9. RAW Spin Recovery Cell Map and Flight Test Trajectory 



Figure 10. Spin Recovery Trajectory Comparison 

It is interesting to compare the RAA and RAW recovery strategies from the perspective of a human pilot in a 
full-scale airplane. The RAA control law ultimately led to a complete recovery, but large excursions in the motion 
variables were encountered before the aircraft recovers. The RAW control law rapidly moves the aircraft from the 
spin critical point to a new critical point on the <7j = 0 switching surface. This strategy is probably preferred by 
pilots because yaw rate is eliminated first, which provides an immediate cue to the pilot that the recovery will be 
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successful. An almost pure rolling motion at zero angle-of-attack defines the new critical point. Once the aircraft 
has reached this new critical point, the pilot completes the recovery using ailerons to stop the rolling motion. This 
second step represents the sliding mode of the VSC control scheme. 

VI. Conclusion 

The value of an integrated approach to modeling, analysis, and control system design lies in the ease at which 
end-to-end flight dynamic studies can be conducted. The nonlinear lifting line modeling method yields a complete 
aerodynamic model yet requires only a small amount of basic geometric information. The subspace cell-to-cell 
mapping method gives a visual picture of the global stability properties near the chosen subspace. The dynamics of 
the state-space are readily manipulated using the theory and methods of variable structure control. The result is an 
integrated computational environment that can be used by a single investigator to rapidly explore flight dynamic 
issues related to large-amplitude maneuvers. 

The end-to-end approach was demonstrated by investigating spin recovery control strategies for a small, hand- 
launched glider. This investigation showed that the recovery path could be manipulated by placing critical points at 
desired locations within the motion state-space. Variable structure control strategies were effectively used to shape 
the flow of trajectories to the critical point. The investigation confirmed that a conventional spin recovery technique, 
which requires the ailerons to be deflected with the spin, is a two-phase process that closely resembles a variable 
structure control law. The reaching mode phase defines a stable switching surface at zero yaw rate. During this 
phase, a simple relay control law drives the aircraft to a critical point that is almost pure rolling motion. The sliding 
mode phase control law would then eliminate the rolling motion, leading to a complete recovery. Flight test 
experiments verified the stability of the reaching mode phase. 
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